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Metal-insulator transition in the two-orbital Hubbard model at fractional band 

fillings: Self-energy functional approach 

Kensuke Inaba 1 and Akihisa KOGA 2 



o 
o 

(N 
C 

Hi 



i 

CO 
-4-* 



I 

C 

o 
o 



> 

00 

m 

o 
l> 
o 



1 Department of Applied Physics, Osaka University, Suita, Osaka 565-0871 , Japan 
2 Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

(Received February 5, 2008) 

We investigate the infinite-dimensional two-orbital Hubbard model at arbitrary band fillings. 
By means of the self-energy functional approach, we discuss the stability of the metallic state 
in the systems with same and different bandwidths. It is found that the Mott insulating phases 
are realized at commensurate band fillings. Furthermore, it is clarified that the orbital selective 
Mott phase with one orbital localized and the other itinerant is stabilized even at fractional 
band fillings in the system with different bandwidths. 
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1. Introduction 

Recently, strongly correlated electron systems with or- 
bital degeneracy have attracted much interest. 1 Among 
them, heavy fermion behavior in the transition metal 
oxides is one of the most important issues in the con- 



£h ' densed matter physics. In the lithium vanadate L1V2O4 



geometrical frustration originating from the pyrochlore 
structure suppresses magnetic correlations, realizing a 
heavy fermion state below T ~ 40K. 3,4 The orbital de- 
generacy in the 3d bands was also suggested to stabi- 
lize the heavy fermion state. 5-7 Other interesting exam- 
ples are the transition metal oxides Ca2- x Sr x Ru04 8 and 
La n+ iNi n 03n+i- 9, 10 It is suggested that a Mott transi- 
tion in some of relevant orbitals, i.e. an orbital selective 
Mott transition (OSMT), 11 is induced by the chemical 
substitution of Ca ions and also by the change in tem- 
peratures, where the heavy metallic state is realized. Fur- 
thermore, the pressure-induced OSMT is also suggested 
in the vanadium oxide V2O3, 12 which stimulates further 
theoretical 13-34 and experimental 35,36 investigations on 
the effect of electron correlations in the system with or- 
bital degeneracy. 

The two-band system in infinite dimensions is one 
of the simplest models with orbital degeneracy, where 
many groups have discussed the nature of the Mott 
transition, 37-49 the magnetism, 50-52 the thermodynamic 
properties, 53,54 etc. Among them, the possibility of the 
OSMT related to the above compounds has been dis- 
cussed recently. 11, 12, 16-18 It has been clarified that sepa- 
rate Mott transitions occur in the half-filled system with 
different bandwidths, in general. 19 On the other hand, 
in a certain parameter region, the effect of the different 
bandwidths is diminished, where a single Mott transition 
is induced. 14,19,23,24 In contrast to these intensive stud- 
ies of the bandwidth control Mott transitions, little is 
known about the filling control Mott transitions, which 
may be important to discuss the stability of the OSMT 
in real materials. Therefore, it is desired to discuss the 
filling control Mott transitions as well as the bandwidth 
control Mott transitions systematically. 

Our concern is to investigate the two-orbital Hub- 



bard model at arbitrary band fillings. For this pur- 
pose, we make use of the self-energy functional approach 
(SFA), 55-57 which allows us to discuss the Mott transi- 
tion systematically. By calculating several physical quan- 
tities, we discuss the stability of the metallic state in the 
system. 

This paper is organized as follows. In Sec. 2, we intro- 
duce the two-orbital Hubbard model. In Sec. 3, we treat 
the two-orbital Hubbard model with same bandwidths 
to test the validity of our analysis. We discuss the na- 
ture of the Mott transition in the system with different 
bandwidths to determine the phase diagram in Sec. 4. A 
brief summary is given in Sec. 5 

2. Model and Method 

We investigate the two-orbital Hubbard model with 
different bandwidths, which is explicitly given as, TL = 
Ho + W, W = £. K with 



Tin — 



(*a-M<*»i)4 
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(1) 
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a. gg' 

- ^(4iT c m4>i c »2T + c\^c\ n c m c ai + H.c.),(2) 

where c\ aiJ (ci aa ) is a creation (annihilation) operator of 
an electron with spin <j(=|,J.) and orbital a(— 1,2) at 
the ith site, and rii aa is the number operator. Here, t a 
denotes the hopping integral for the ath orbital, /i the 
chemical potential, U (If) the intra-orbital (inter-orbital) 
Coulomb interaction, and J the Hund coupling including 
the spin- flip and pair- hopping terms. We impose the con- 
dition U = U' + 2 J, which is obtained by the symmetry 
arguments for the degenerate orbitals. 

To discuss the competition between the metallic and 
the Mott insulating phases, we make use of the SFA pro- 
posed recently. 55, 56 Since this method is based on the 
variational principle, it has an advantage in discussing 
the nature of the Mott transitions. In fact, it has been ap- 
plied to correlated electron systems at half filling, where 
the precise phase diagrams have been obtained. 30, 54-57 
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Here, we deal with the two-orbital system at arbitrary 
band fillings to determine the phase diagrams. The de- 
tail of the SFA for the doped system is explicitly shown 
in Appendix. In the paper, we use a semi-circular density 
of states (DOS), p a (u) = 4/irW a y/1 - (2u/W a ) 2 , where 
W a is a bandwidth for the ath orbital. 

In the next section, we treat the degenerate Hubbard 
model with same bandwidths as a simple model, to check 
the validity of our analysis for the doped system. Then 
we determine the phase diagram in the system to discuss 
the role of the Hund coupling for the Mott transitions at 
quarter filling. 

3. Mott transitions in the system with same 
bandwidths 

We consider the two-orbital Hubbard model with same 
bandwidths (W\ = W2 = 4). The renormalization fac- 
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Fig. 1. Renormalization factor as a function of the total electron 
count n in the case J/U = (a) and J/U = 0.1 (b). 



tor Z, which is proportional to the inverse of the ef- 
fective mass, is shown in Fig. 1. First, we focus on the 
case J/U = 0, as shown in Fig. 1 (a). The introduction 
of the Coulomb interaction monotonically decreases the 
renormalization factor and the double dip structure ap- 
pears in the case (U > 4). In the strong coupling region, 
the renormalization factor reaches zero when n — 1 and 
n = 2, while it never vanishes otherwise. 40,41 This result 
is consistent with the well-known fact that the Mott tran- 
sition occurs only at commensurate band fillings such as 
quarter filling (n = 1) and half filling (n = 2). 40, 41 - 48 
The corresponding critical points U c <~ 9.2 and 7.7 are in 
agreement with those obtained by other numerical meth- 
ods. 45 ' 48 Therefore, we can say that the reliable results 
are obtained in terms of the SFA not only at half filling 
but also at arbitrary band fillings. A similar dip struc- 
ture appears in the system with a finite Hund coupling 
J/U = 0.1, as shown in Fig. 1 (b). 

We now address how the Hund coupling affects the 
competition between the metallic and the Mott insu- 
lating phases at quarter filling (n = 1), in comparison 
with the results for the half filled case (n — 2). 48 ' 49 ' 54 
In Fig. 2, we show the renormalization factor, and the 



spin and orbital susceptibilities. In the case, the param- 
eter U = U' — J is useful to discuss the Mott transi- 
tion in the system, which corresponds to the energy scale 
of the lowest Hubbard gap [see the inset of Fig. 2 (a)[. 
When J/U is small, the increase of the interaction U 
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Fig. 2. The quasi-particle weight (a), the spin susceptibility (b), 
the orbital susceptibility (c) as a function of U when n = 1. The 
inset of the panel (a) is the density of state — lmG a (w)/n for 
J/U = 0.2, U = 15 and U' = 12, where both the renormalizcd 
metallic state at Fermi level and the precursor of the Hubbard 
satellites appear. 



decreases the renormalization factor, where spin and or- 
bital susceptibilities are monotonically enhanced. On the 
other hand, the large Hund coupling (J/U = 0.3) leads 
to different behavior. It is found that the introduction 
of U rapidly decreases the renormalization factor, where 
the spin (orbital) susceptibility is increased (decreased). 
When U = U' - J < W <C U,U' (U = 0.2), the Hub- 
bard gap does not appear near the Fermi level, but high- 
energy structures appear due to U and U', as shown in 
the top panel of Fig. 3. It suggests that the spin-triplet 
two-electron state formed by J plays an important role in 
low energy states, while the spin-singlet states contribute 
the high energy features. However, when the system ap- 
proaches the Mott critical point (the bottom panel of 
Fig. 3), the Hubbard satellite that originates from the 
spin-triplet state appears at the position of U, while the 
Hubbard satellites have roots in the spin-singlet states 
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Fig. 3. The density of states for J/U = 0.3 with several choices 
of U = 0.2, 2, and 5. 



are shifted to higher energy region. Then, the triplet 
state does not directly contribute to low energy states, 
but contributes through a virtual process between the 
ground state and the triplet state. Therefore, the heavy 
fermion state is gradually changed in character around 
U ~ 0.3. Eventually, both susceptibilities simultaneously 
diverge at a critical point U c — 5.9 (U c = 58.7), where 
the second-order transition occurs to the Mott insulating 
phase, as shown in Fig. 2. 

By performing similar calculations, we end up with the 
phase diagram as shown in Fig. 4. Note that at quarter 
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Fig. 4. The phase diagram for the system with same bandwidths 
at quarter filling. 



filling, spin and orbital fluctuations are simultaneously 
enhanced at the critical point, yielding the second-order 
Mott transition. This behavior is in contrast to that for 
the half filled case. 45, 48,54 At half filling, the second-order 
Mott transition occurs in the case U = 0, where spin and 
orbital fluctuations are enhanced simultaneously. On the 
other hand, if the Hund coupling is introduced, orbital 
fluctuations are strongly suppressed, yielding the first- 
order Mott transition. 32, 37,48 

In this section, we have studied the degenerate Hub- 



bard model with same bandwidths by means of the SFA. 
The obtained results reproduce the well-known numer- 
ical results, 40,41,48 implying that the SFA allows us to 
investigate strongly correlated electron systems with and 
without particle-hole symmetry systematically. In the 
following, we consider the two-orbital system with dif- 
ferent bandwidths to discuss ground state properties at 
arbitrary band fillings. 

4. The system with different bandwidths 

Let us move our attention to the effect of different 
bandwidths. Recently, theoretical advances have been 
made in the system at half filling, where the nature of 
the OSMT has been discussed. 14,19,23 - 26,28,31,34 It has 
been clarified that separate Mott transitions occur in 
general, which merge to a single Mott transition under 
a certain condition. 14,19,23,24 Furthermore, it has been 
discussed how stable the OSM phase is against the hole 
doping, 19 the hybridization, 20 the magnitude and/or the 
anisotropy of the Hund coupling, 15, 19,30 etc. However, it 
is not trivial whether the OSM phase in the doped system 
is adiabatically connected to the Mott insulating phase 
at quarter filling. This may be important to understand 
real materials, in which the total electron count can be 
tuned by the chemical substitutions. Therefore, it is nec- 
essary to systematically discuss how the Mott insulating 
phases, which are realized in commensurate band fillings, 
are stabilized in the doped system, which may be a key 
to understand the Mott transitions in real materials with 
orbital degeneracy. 

Here, we consider the Mott transitions in the Hub- 
bard model with the different bandwidths W\ = 2 and 
W2 = 4 (R = 0.5). To overview the competition between 
the metallic and the Mott insulating phases, we first cal- 
culate the renormalization factor for each orbital in the 
system with J/U = 0.1, as shown in Fig. 5. When n = 2, 
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Fig. 5. Renormalization factor for narrower (wider) orbital (a) 
[(b)] as a function of the total electron count n in the system 
with R = 0.5 and J/U = 0.1 



half filling is realized in each band, where the effective 
Coulomb interaction should depend on the bandwidth 
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of each orbital (U%ff — U/W a ). Therefore, with the in- 
crease of the interaction, the renormalization factor for 
the narrower band are decreased rapidly. Eventually, the 
OSMT occurs at the critical point U c <~ 3.3, inducing 
the OSM phase with one orbital localized and the other 
itinerant. Further increase of interaction induces another 
Mott transition in the wider band at the critical point 
U c ~ 4.0. These results are consistent with the previous 
results. 24 ' 30 In contrast, somewhat complicated behavior 
appears away from half filling, where an electron count 
for each band depends on the magnitude of the interac- 
tions. 

To make this clear, we calculate the renormalization 
factor and the electron count for each orbital with a 
fixed total electron count n = 1.0 (quarter filling), 1.1 
and n = 1.6, as shown in Figs. 6 (a), (b), and (c). Now, 
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Fig. 6. The renormalization factor for narrower (wider) orbital 
(a) [(b)] as a function of the Coulomb interaction U in the system 
with R = 0.5 and J/U = 0.1 



0.8 



0.6 



0.4 



0.2 



(a)w=l 



(b) n=i.l 



(c)n=1.6 




2 4 6 
U 



20 40 60 2 4 6 8 1012 

U U 



Fig. 7. The electron count for narrower (wider) orbital (a) [(b)] 
as a function of the Coulomb interaction U in the system with 
R = 0.5 and J/U = 0.1 



we focus on the case n = 1.6, where the total electron 



count is close to half filling. As increasing the interac- 
tions, the renormalization factor for the narrower band 
reaches zero at a critical point U c = 11.7, where half fil- 
ing is realized in the narrower band. This implies that 
the commensurability for the narrower band is satisfied 
in the presence of electron correlations, and thus the sys- 
tem is driven to the OSM phase. 19 On the other hand, 
when the system is close to quarter filling (n = 1.1), the 
electron count for both bands is far from unity even in 
the strong coupling regime, as shown in Fig. 7 (b). This 
implies that the heavy quasi-particle state persists and 
the Mott transition never occurs (at least, up to U <~ 60). 

In contrast, at quarter filling (n = 1), the renormaliza- 
tion factors z\ and z-i reach zero simultaneously, where 
the single transition occurs. It is found that the electron 
count for each orbital still remains fractional even in the 
Mott insulating state. This suggests that the commensu- 
rability is never satisfied for cither orbital. An important 
point is that the Mott transition at quarter filling is quite 
different from the OSMT discussed above, implying that 
the Mott insulating state is not adiabatically connected 
to the OSM state. To clarify the difference in these Mott 
transitions, it is instructive to discuss spin and orbital 
fluctuations in the system at arbitrary band fillings. 

In Fig. 8, we show the local susceptibilities for spin 
and orbital degrees of freedom. Characteristic behavior 
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Fig. 8. Spin (a) and orbital (b) susceptibilities as a function of 
the total electron count n when J/U = 0.1. 



of the Mott transition appears in the figure. When n <~ 1, 
the increase of the interaction enhances spin and orbital 
fluctuations simultaneously, inducing the Mott transition 
at quarter filling. On the other hand, when the system 
is close to half filling (n > 1.3), the spin susceptibility is 
enhanced and the orbital susceptibility is strongly sup- 
pressed in the strong coupling region, where the OSM 
phase is realized with one orbital localized and the other 
itinerant. An important point is that the metallic state 
between these Mott insulating phases are quite sensitive 
to the total electron count, as discussed above. Therefore, 
it is expected that the ground state is drastically changed 
by some perturbations, e.g. the crystalline electric field, 
the chemical substitution, etc. 
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We end up with the phase diagram as shown in Fig. 9. 
At half filling (n = 2), as introducing the Coulomb inter- 
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Fig. 9. Phase diagram for the system with different bandwidths 
R = 0.5 when J/U = 0.1 (main) and J/U = 0.25 (inset). The 
shaded region represents the OSM phase with one orbital local- 
ized and the other itinerant. 

action, the Mott transition in the narrower band occurs 
to the OSM phase. Further increase of the interaction 
drives the system to the Mott insulating phase. If the 
total electron count is changed in the Mott insulating 
phase, holes are doped only in the wider band due to 
the strong correlations in the narrower band. Therefore, 
the OSM phase is stabilized even away from half filling. 
When the electron count approaches quarter filling, the 
holes are introduced in both bands, thus the system is 
driven to the metallic state. In contrast, the different 
type of the Mott insulating state appears beyond a cer- 
tain critical point when the system is quarter-filled. 

Up to now, we have restricted our discussion to the 
condition J/U = 0.1 and R = 0.5. If the Hund coupling 
is increased, orbital fluctuations are further suppressed. 
Therefore, it is expected that the OSM phase becomes 
more stable, while the Mott insulating phase at quar- 
ter filling becomes unstable. In fact, in comparison with 
the phase diagram for J/U = 0.1 and R = 0.5, the re- 
gion of the OSM phase broadens and the critical point 
at quarter filling is increased, as shown in the inset of 
Fig. 9. The large difference of the bandwidths also sup- 
presses orbital fluctuations, as discussed in the previous 
paper. 32 Therefore, the same phase diagram is expected 
to be obtained, where the wide OSM phase away from 
half filling is realized. 

5. Summary 

We have investigated the two-orbital Hubbard model 
with same and different bandwidths at zero temperature. 
By making use of the SFA, we have discussed the com- 
petition between the metallic and the insulating states 
in the system at arbitrary band fillings. We have found 
that the Mott transition occurs only at commensurate 
band fillings in the system with same bandwidths. On the 
other hand, in the system with different bandwidths, the 
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OSM phase with one orbital localized and the other itin- 
erant is stabilized in a certain parameter region. There 
are some open problems in the present study. One of 
the most important issues is magnetism of the system, 
which has not been treated here, since we have restricted 
our attention to the paramagnetic phase. This problem 
is under consideration. 
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Appendix: Formulation of the SFA 

We briefly summarize the formulation of the SFA in or- 
der to discuss the Mott transition in the system without 
particle-hole symmetry. Here, by combining the numer- 
ical method with the analytic method, we describe the 
Mott insulating state in the framework of the SFA. 

First, let us consider the free energy F for the system 
as, 

F/L = Q/L + I m, (A-l) 

where L is the number of sites and n is the electron count 
per site. The grand potential O is given as, 

fi[£] = r[S] + Trln[-(G- 1 -S)- 1 ], (A-2) 

where T[S] is the Legendre transformation of the 
Luttinger-Ward potential. 62 G and £ are the bare 
Green function and the self-energy, respectively. The 
condition dfl[S]/d'E = gives us the Dyson equation 
G 1 = G^ 1 - £, 62 where G is the full Green func- 
tion. Note that the potential does not depend on 
the detail of the non-interacting Hamiltonian. 55 This en- 
ables us to introduce a reference system with the same 
interacting term. The Hamiltonian is explicitly given by 
W re f(t') = Ho(t' ) + H' with the parameter matrix t'. 
Then the grand potential is given as, 

n[s(t')] = fi(t') 

+ Trln[-(w + ^-t-S(t')) _1 ] 

- Trln[-(w + /i-t'-S(t')) _1 ](A-3) 

where O(t') and S(t') arc the grand potential and 
the self-energy for the reference system. The condition 
dfl[E(t')]/dt' = gives us an appropriate reference sys- 
tem in the framework of the SFA, which approximately 
describes the original correlated system. 

Here, we introduce a two-site Anderson impurity 
model as the reference system, which has successfully 
described strongly correlated electron systems in the in- 
finite dimensions. 30 ' 54 " 56 ' 63 - 64 The Hamiltonian of the 
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reference system is explicitly given as 



H 



rcf — 



(A-4) 



H 



W _ 



rcf 



H.c.)+n' t , (A-5) 



where aJ Q(T (aic(CT) creates (annihilates) an electron with 
spin cr and orbital a, which is connected to the ith site 
in the original lattice. The grand potential per site is 
rewritten as, 



+ 



2 EE 



^Pa(%km(z)], (A-6) 



5(2;) = -Tln[l + exp(-a;/T)] 



(A-7) 



where f2i mp is the grand potential for the reference sys- 
tem and Lu am (z) [oj' am ] are the poles of Green function 
G (G 1 ) for the original (reference) system. The Green 
function of the reference system is given as, 

G» - [u> + M ~ e% - J»_ ra - S„M](A-,8) 

UJ -\- p £ a 

where S a (w) is the self-energy for the ath orbital. On 
the other hand, the full Green function for the original 
system is given as, 

G a (w) = j dzp a (z)G a (uj;z), (A-9) 

G a {u;z) = [w + n-z-Haiu)]- 1 , (A-10) 

where p a {x) is the density of states for the bare ath 
orbital. Note that the following differential equation is 
efficient to deduce the poles uj am (z) of the Green function 
for the original system, 



du) am (z) 
dz 



duo 



(All) 



By solving these differential equations, we estimate the 
free energy [Eq. (AT)] numerically. The variational pa- 
rameters A = (p, V a , e£, e£) are determined by the con- 
dition dF/dX = 0, which corresponds to the stationary 
point in the parameter space. 

In the framework of SFA, the variational parameters 
for the metallic and the OSM states can be determined 
numerically. In contrast, it may be difficult to determine 
unique variational parameters for the Mott insulating 
state, as discussed above. To overcome this, we make 
use of the analytical and the numerical methods. 

For simplicity, we consider the two-orbital model with 
J/U = at quarter filling. Recall that the Mott insulat- 
ing phase is simply represented in the atomic limit, by 
taking into account local electron correlations. Since the 
four one-electron states c4 CT |0) are degenerate (ef = ejj), 
the ground state |$g) m the reference system [eq. (A-5)] 
should be given as 



|4> G >=E 



10), 



(A12) 



where n a is the expectation value of band fillings for the 
ath orbital, which satisfies the constraint J2 aa . n a = 1. 

The Green function and the self-energy for the refer- 
ence system [eq. (A-5)] are given as 

1 - n a 



E a (w) = 



+ 



--(A 13) 



uj + p + e c a uj + fi + s c a - U 

uj + p + el-l/G'Juj). (A- 14) 

The Green function for the original system is then writ- 
ten as G a (uj) = J dzp a (z)G a (uj; z) with 

G a (uj;z) = 



uj + p — z — S a (w) 



E 

i=± 



(A-15) 



where 



z + s c a + U(l-2n a ) 



RUz) 



1± 



r](z) 



V(z) = VU 2 + 2(1 - 2n a )(z + e c a )U + (z + e%f . 

This formula is equivalent to the Hubbard-I type approx- 
imation. Finally, we obtain the free energy per site for 
the system as 

^/ L = E^ +e «)( 1 " n «)+^E n « 

a a 

-EE / dz Pa {z)g{z)0[-tf{z)]. (A-16) 



+ 



a l=± • 



At the stationary point in the free energy, three con- 
ditions are obtained as, 

1 = E / dzp a {z)6[-^{z% (A- 17) 

i=± J 

n a = E / dz Pa {z)R?{z)6[-ti?{z)l (A- 18) 
i=± J 

We wish to note that the right hand side of eq. (A- 18) 
is equivalent to G a (uj)dbj/2Tri for each orbital. This 
implies that n a is identical to the electron count for the 
lattice system. 

As described above, the chemical potential is difficult 
to be determined uniquely in the insulating phase. In 
fact, the condition dF/dp = is satisfied in the case 
WJ2 < p < U - W a /2 {U » W a ), which originates 
from the presence of the Hubbard gap. Here, we set the 
chemical potential as p — U/2. Then the condition eq. 
(A- 17) is always satisfied, and eq. (A- 18) is rewritten as 



n a = J dzp a 



(z)R a _(z). 



(A-19) 



By solving the self-consistent equations eq. (A-19), 
^2 aa n a = 1, and e\ — e\, we can describe the Mott 
insulating phase in the framework of the SFA. 
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